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ABSTRACT 

We have measured velocity dispersions (a) for a sample of 36 galaxies with J < 21.2 or M r < —20.6 
mag in MS 1054-03, a massive cluster of galaxies at z — 0.83. Our data are of uniformly high quality 
down to our selection limit, our 16-hr exposures typically yielding errors of only <5(<r) ~ 10% for L* and 
fainter galaxies. By combining our measurements with data from the literature, we have 53 cluster 
galaxies with measured dispersions, and -ffST/ACS-derived sizes, colors and surface brightnesses. This 
sample is complete for the typical L* galaxy at z ~ 1, unlike most previous z ~ 1 cluster samples which 
are complete only for the massive cluster members (> 10 11 M@). We find no evidence for a change in 
the tilt of the fundamental plane (FP). Nor do we find evidence for evolution in the slope of the color-cr 
relation and M/Lg-a relations; measuring evolution at a fixed a should minimize the impact of size 
evolution found in other work. The M/Lb at fixed a evolves by Alog 10 M /Lb = —0.50 ±0.03 between 
z = 0.83 and z = 0.02 or d log 10 M/L B = -0.60 ±0.04 dz, and we find A(U-V) Z = -0.24 ±0.02 mag 
at fixed a in the rest-frame, matching the expected evolution in M/ Lb within 2.25 standard deviations. 
The implied formation redshift from both the color and M/Lb evolution is = 2.0 ± 0.2 ± 0.3(sys), 
during the epoch in which the cosmic star-formation activity peaked, with the systematic uncertainty 
showing the dependence of z+ on the assumptions we make about the stellar populations. The lack of 
evolution in either the tilt of the FP or in the M/L- and color-cr relations imply that the formation 
epoch depends weakly on mass, ranging from = 2.3IJ3 at cr = 300 km s _1 to = 1.7^q 2 at 
a = 160 km s _1 and implies that the inital mass function similarly varies slowly with galaxy mass. 
Subject headings: galaxies: clusters: general — galaxies: elliptical and lenticular, cD, — galaxies: 

evolution — galaxies: fundamental parameters — galaxies: photometry — clusters: 

individual: EMSS 1054.4-0321 



1. INTRODUCTION 

The evolution of the mass-to-light ratio (M/L) of 
early-type galaxies can be measured by combining the 
velocity dispersion (cr), the effective radius (r e ), the av- 
erage surface brightness within the effective radius ((I e )), 
with which we can measure M/L cx a 2 /(r e (I e )) and 
how it depends on mass, M cx r e a 2 . Typically, these 
three variables are combine d into an empirical re lation 
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The absolute M/L value and the how fast the M/L 
evolves with time can be used as a technique for mea- 
suring the age of stellar population (|Tinsfevl [r972'). The 
younger the population, the more luminous the stars, 
and so it will have lower M/L values and the M/L will 
evolve faster. Massive cluster galaxies, those galaxies 
with M > 10 11 M Q where M cx cr 2 rv., out to z > 1 ap- 
pear to evolve as A In M/L b — z (Ivan Dokku m fc Franx 
1996 1 : iKelson et all [19971 l2000d: Ivan Dokkum fe Franx 
200ll Ivan Dqkkum fe Stanford 120031; iWuvts et all 120041: 
Holden et all l2005bl: Ivan Dokkum fc van der Marell 
2007I K This rate of evolution 
field sample s, though there is 
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1 20031: Ivan Dokkum fc Ellisl 120031: Ivan de Ven et al 
pOOa Ivan der Wel et all 120041: iTreu et al.l l2005allbl : 
Ivan der Wel et aLll2005f) . This rate of evolution implies 
an epoch of formation for the stars in early-type galaxies 
of z+ ~ 2, assuming a passively evolving simple stellar 
population with a standard Salpeter-like initial mass 
function (IMF). This redshift is the near the peak of star 
formation ( e.g; iMadau et al.lll998l: [Steidel et al.lH99l 
iDickinson etall 120031: iRudnick et al.l 120031 ) . Thus, the 
majority of stars in massive early-type galaxies formed 
around the time that the average star in the universe 
was formed. 

The FP results of lTreu et all (|2005bL hereafter T05), 
find that for galaxies with masses of < 10 11 Mq, the 
M/L values are lower and evolve more rapidly than 
those for higher mass galaxies. This implies a typi- 
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cal luminosity weighted age of the stellar populations 
co rresponding to a redshift between z+ = 1 and = 
2. Ivan Dokkum fc van der Marell (|2007l hereafter vv07) 
confirmed that at redshifts of z ~ 1, the FP of clus- 
ter galaxies shows a different distribution of M/L val- 
ues than at z ~ with some lower mass galaxies 
(< 10 11 Mq) having low M/L and potentially more re- 
cent epochs of f o rmation than the higher mass galaxies. 
|j0rgensen et al.l (|2006l l2007f ) found a correspondingly 
steeper tilt to the FP for cluster galaxies in RX JO 15 2- 
13, indicating that cluster galaxies show similar behav- 
ior to field galaxies, though with a smalle r magnitude 
effect. Similar results are clearly seen in iSaglia et al.l 
(2010) which jointly studies a field and cluster sample 
with a uniform selection. The observed increase in the 
tilt of the FP was interpreted as a younger average stel- 
lar population age. However, for a sample limited in L, 
only tho se galaxies with a l ow M /L will appear at low 
masses. Ivan der Wei et all (|2005l vdW05) showed that 
it is difficult to distinguish between evolution in the tilt 
and in the scatter of the FP. This is the direct result of 
the selection by optical magnitude in their sample. They 
conclude that the z ~ 1 FP is likely to be different from 
the present-day FP, but that the form (tilt or scatter) of 
the evolution is not constrained besides the evolution in 
the zero-point at high mass. 

Younger stellar populations will not only have a lower 
M/L, but also have bluer colors. As we look to higher 
redshift stellar populations where the slope of the FP ap- 
pears to tilt, we should see a corresponding trend of bluer 
galaxies at lower mass, or a change in both the slope of 
the color-magnitude and color-mass relations, with the 
cavea t that magnitud e and m ass are not the same quan- 
tity. iBlakeslee et~aTl (|2006l hereafter B06) found that 
the mean and scatter in the color-magnitude relation of 
cluster early-type galaxies at z = 0.83 corresponded to a 
formation epoch of z* = 2.2 ove r a broad r ange o f mag- 
nitudes, and therefore masses. IMei et al.l ((2001 found 
no evolution in the slope of the color-magnitude relation 
in cluster galaxies out to z = 1.3 extending the trend 
found with lower r edshift clusters ([Stanford et al.lll998t 
IHolden et al.l 120041). Th is is complementary to th e re- 
sults of [Bell et al.l pOOl and lRuhland et al.l(|2009t) who 
found little evidence for evolution of the slope of the 
color-magnitude relation of field galaxies. A lack of evo- 
lution in the slope of the color-magnitude relation when, 
at the same time, lower mass galaxies are to apparently 
h ave lower mea n M / L b values is puzzling. 

Ivan Dokkuml (|2008l . hereafter vD08) combined the 
M/ L evolution and the evolution in the U — V color of 
cluster early-types to measure both the epoch of galaxy 
formation for massive galaxies and to constrain the slope 
of the initial mass function (IMF). Evolution in the col- 
ors of passively evolving stellar populations are governed 
by the location of the main-sequence turn-off, which is 
determined by the age of the stellar population. The rate 
of M/L evolution, however, is determined by the slope 
of the IMF at the turn-off as well a s how bright t he stars 
at the main-sequence turn-off are ((Tinslcy 1972). Thus 
a slow pace of color evolution and a fast pace of M/L 
evolution, for example, can be explained by modifying 
the initial mass function. vD08 found, for the more mas- 
sive cluster galaxies (> 10 11 M@), a rapid pace of M/L 
evolution as compared to the color evolution. Thus, he 



concluded that the slope of the initial mass function at 
the main-sequence turn-off was flat, x = —0.1 instead of 
the usual Salpeter x = 1.35. This surprising result did 
not match previous measurements over the same red- 
shift range made by other methods (Kelson et al. 2000c, 
I200in . a difference that c ould be explained by sample 
selection or the improved Marastonj (|2005f ) population 
synthesis models used by vD08. This result is also in 
conflict with some local indirect measurements of the 
IMF through lensing (eg. iTreu et al.l 120101: lAuger et all 
I2010T) . though possible accommodation can still be made 
through variation in the dark matter halos properties. 
Because of the IMF, vD08 finds that the stars formed 
with a earlier formation epoch, z — 3.7lg'|. vD08 pro- 
poses a model that, by z — 2, predicts a more normal 
IMF. The rapid M/L evolution in the FP of lower mass 
galaxies implies that they should have formed at about 
z = 1 — 2. However, the lack of evolution in the slope of 
the color-magnitude relation would imply an even steeper 
IMF slope for low mass galaxies than vD08 found for 
the high mass galaxies. This is in marked contrast to 
the expectations from the mode l of vD08 or othe r IMF 
slope measurements near z — 2 (jBlain et al.lfl999h . both 
of which expect a more normal IMF for galaxies with 
masses < 10 11 M Q . 

Because of the magnitude limits used for selection in 
previous work, computing the typical stellar age of for- 
mation for < 10 11 Mq galaxies has been challenging. The 
highest M/L galaxies, those with the older populations, 
cannot be directly observed at < 10 11 M Q and the im- 
pact of these galaxies must instead be included through 
a modeling process as was done in T05 or vdW05. We 
selected a sample of z ~ 1 early-type galaxies in the field 
of MS 1054-03, a cluster of galaxies at z = 0.83 that sam- 
ples much farther down the mass function. The main goal 
was to search for the more rapid evolution of the lower 
mass galaxies implied by the results of other efforts. A 
summary of the sample selection, and our measurements 
of the parameters for the FP, the effective radius (r e ), 
the average surface brightness within the effective radius 
((/ e )), and the velocity dispersion (a) are discussed in 
Section [2l We determine the completeness of our sam- 
ple in Section [3l Using these measurements, appropri- 
ately weighted, we measure the FP, color evolution, and 
the Mj Lb evolution in Section [5] We use these mea- 
sures to compute the typical star-formation epoch of the 
stellar populations and comment on constraints of the 
IMF, in Section [5l Our results are summarized in Sec- 
tion [6l Throughout this paper, we assume £l m = 0.27, 
Q A = 0.73 and H a = 71 km s" 1 Mpc" 1 . 

2. DATA AND MEASUREMENTS 

We have a sample of 59 galaxies in the field of the 
z = 0.83 cluster MS 1054-03 for which we have obtained 
spectra. Of these, 36 have dispersion measurements. 
For these galaxies we have high resolution imaging using 
the Advanced Camera for Surveys (ACS) on the Hubble 
Space Telescope (HST). Using the HST/ ACS images, we 
measure the sizes and surface brightness by fitting ellip- 
tical Sersic models to the data. The velocity dispersions 
are measured by fitting stellar templates broadened by 
the appropriate Gaussian line of sight velocity disper- 
sions, representing the typical motion of stars inside the 
galaxies. Below we discuss in detail our observations and 
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Figure 1. Example spectra and images. All images are 5'.'0 on a side and all spectra are in the observed frame. The vertical gray band 
masks the atmospheric absorption from the A band. The top three galaxies represent typical objects in our sample, 19.5 < J < 20.5 with E 
and SO morphologies and logo - > 2.2. The brightest of these galaxies, w5756, is one of the three galaxies from W04 which we re-observed. 
5234 is the faintest galaxy for which we measured a a value. The bottom two galaxies are both late-type systems, with the stretch for 
7534 adjusted to illustrate this. 7534 has strong enough metal absorption lines for us to measure a dispersion. In 11461 the strong Balmer 
absorption and lack of metal lines prevented us from measuring a cr, despite its bright J magnitude. Both 5234 and 11461 were imaged 
using the /814 on ACS while the remaining galaxies were observed with 1775 filter. 



the resulting measurements. 

2.1. z — 0.83 Sample Definition 

We based our sample on a magnitude limited survey 
of galaxies in the field of MS 1054-03. The original se- 
lection was done using a limiting magnitude of J = 21.2 
mag (Vega) from a catalog of galaxies observed with the 
CFH12K at CFHT and the WIRC camera on the du 
Pont. Each galaxy in this sample had either a spectro- 
scopic redshift, or a photometric redshift from a catalog 
of BVRIJH imaging data. Galaxies with spectroscopic 
redshifts were preferred in our selection. These photo- 
metric redshifts, the catalog construction, and the spec- 



troscopic subsample will be discussed in a later paper. 

Beyond meeting the J magnitude limit, the main crite- 
rion for selection was that the galaxy lie in the ACS field 
of view. We selected a DEIMOS mask center and posi- 
tion angle that maximized the number of galaxies with 
ACS imaging and with J < 21.2. A total of 44 galaxies 
met these criteria. Beyond our main sample, 12 supple- 
mental galaxies were selected to have 21.2 < J < 21.8 
and an additional four were targeted outside of the HST 
field of view. Also, we re-observed thr ee galaxies that 
were m the sample of lWuvts et all (|2004 hereafter W04). 
Therefore, the total number of galaxies targeted for dis- 
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Figure 2. Color-magnitude relation for galaxies in the z = 0.83 cluster MS 1054-03 using HST/ ACS data. We plot the Veoe— *775 color 
of galaxies as a function 1775 magnitude, note different y-axis scale for top and bottom panel. Gal axies for which we have Igi4 data, we 
convert the Igi4 and V606 - ^814 data into the 1775 and Veoe— 1775 using the relations in Section |3. II Larger red. filled circles and triangles 
are E and SO galaxies with dispersions while smaller orange diamonds arc E and SO galaxies without dispersions. Green spirals are spiral or 
irregular galaxies, with filled spirals representing those with measured dispersions. We show, with solid lines, the color-magnitude relation 
of B06 while the gray regions shows the 3cr scatter. We will consider all galaxies within that shaded region as on the red-sequence, those 
galaxies with dispersions that lie outside of that region are shown with filled triangles. We find a tight color-magnitude relation for the 
cluster E and SO population in good agreement with the results from B06. This relation extends out to the E and SO galaxies beyond the 
core of the cluster. We note that the statistical errors on the colors are smaller than the symbols used for plotting, see also B06. 



persions was 62. 

We also combine our sample with that from W04. This 
earlier sample has a br ighter I selection , see W04 for de- 
tails, and used LRIS (jOke et al.l 119951 ). with a smaller 
field of view than DEIMOS, for spectroscopic observa- 
tions. The result is that the combined catalog of our 
sample and that of W04 will contain more galaxies in 
the inner regions of the cluster and will have a higher 
completeness at brighter magnitudes. 

Our sample of galaxies is listed in Table [3] in Appendix 
Al. We list our observed and rest-frame quantities along 
with the velocity dispersions. This table includes the 
sample of W04, which are identified by a "w" in front of 
the identification number. 

2.2. HST /A CS Imaging 

The ACS imaging used in this paper comes from two 
separate programs. The first, discussed in B06, used the 
F606W, F775W and F850LP filters to image the center 



of the cluster. We will refer to the filters as Vao6, £775, 
and ^850 for the rest of the paper. The second program 
imaged the outer regions of MS 1054-03 with the V606 
and the F814W, or Igi4 filter. Example images can be 
seen in Appendix A, Figure [1] 

We fit a model to each galaxy to determine the effective 
radius r e and surface brightness (I e ) in the filter that 
most closely matches the rest-frame B. For the inner 
regions, this is the ^775 filter while for those galaxies at 
larger radii, we used Isu- For each galaxy we fit a model 
with a free Sersic paramet er constrained to fall within 
1 < n < 4 using GALFIT ()Peng et all 120021 ) . For this 
paper, we use the best-fitting Sersic parameters for r e 
and (I e ), as this w as shown to y ield the same resulting 
FP as fixing n = 4 (jKelson et alJl2000af ). When the best 
fitting Sersic n was at the limit of the allowed range, i.e. 
n = 1 or n = 4, we refit the image fixing n to that value. 

The total magnitude we use is the normalization of the 
model fit. To measure colors, we used the circularized 
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half-light or effective radius, r e = ahir\fq, where q is the 
ratio of the minor to major axis, or 1 — e, and a^ir is the 
half-light radius along the major axis of the best fitting 
elliptical model as determined by GALFIT. This is the 
radius used both in our M/L analysis and the radius of 
the aperture for which we measured the galaxy's color, 
sec B06 for details. From this circularized r e , we compute 
the average surface brightness with in that radius, or (I e ) . 

Using the simulation results of iHolden et al.l ([2009), 
which placed real galaxy images in the ACS frames, we 
find that our measurements of the total magnitudes are 
too bright by —0.08 mag in the 1775 data or —0.04 mag 
in the I$i4 data. We find a scatter of a = 0.07 mags 
around the input magnitude for the simulations, after 
removing the offset, in good agreement with the expecta- 
tions from total magnitu de measurements in ACS imag- 
ing (jHolden et al.l l2005al ). We apply this offset to our 
data in later sections. 

In addition to fitting a bounded Sersic model, we fit 
a fixed de Vaucouleur's model to our data so we can 
compare the two resulting fits. As expected, for galax- 
ies classified as E or SO the median best fitting n value 
was n = 3.6 while, for those galaxies classified as spirals 
the median n was n = 1.04 (we discuss the classifica- 
tions in the next section.) The n values showed a larger 
scatter for the spirals than for E and SO systems as well. 
In cases where the best fitting n > 3, the values for r e 
are good agreement with the n = 4 values, with differ- 
ences of 11±6% on average. When best fitting bounded 
Sersic n < 3, the differences in the r e are much larger, as 
expected (|Kelson et al.ll2000al) . 

2.2.1. Morphologies 

Galaxy morphologies were obtained from the literature 
for all the galaxies in our sample. All of the morphologies 
were determined by Marc Postman from a J selected cat- 
alog with a limit of J = 22.5 mag. Each galaxy was clas- 
sified in the passband closest to the rest- frame B, either 
the 1775 or the isi4- Th e galaxies were classifi ed in the 
same manner as those in [Postman et al. (2005;|), and the 
galaxies in the central regions of MS 1054-03 were clas- 
sified twice, on c e for this paper and once for the original 
iPostman et al.l ()2005D study. The scatter for the fraction 
of E+S0 galaxies was 6%, the same as was fo und by using 
multiple classifiers in IPostman et al.l (120051) . The mor - 
phologies, using the T system of IPostman etafl ([20"05h . 
are listed in Table [3] Galaxies with T < are E and SO 
systems. 

2.2.2. Redshifted Magnitudes 

We redshift the U, B and V passbands to z=0.83 and 
use the observed photometry to compute what magni- 
tudes we would observe in these redshifted passbands. 
Hence, these passbands measure the flux densities at the 
redshifted effective A for each filter; Aj/(1 + z), A#(l + z) 
and Ay(l + z). For brevity, we will call these redshifted 
magnitudes and label these magnitudes U z , B z and V z . 
With the removal of the additional dimming caused by 
the distance to the galaxies, these magnitudes would be- 
come rest-frame magnitudes. We compute these trans- 
formations between observed a nd redshifted magni tudes 
using a similar proc ess as B06, IHolden et all ((20061 ) and 
IHolden et al.l (|2007l) . but based on the same basic ap- 



proach as Ivan Dokkum fc Franxl (|1996l ) and lKelson et al.l 
(|2000aD . 

Both the 1775 and igi4 filters are mapped into B z . We 
comput ed magnitudes at b o th z = and z — 0.83 for 
a set of IBruzual fc Charlotl ([200l BC03) r model tem- 
plates. For the z = fil ters, we used the U, B and V 
curves from|Buser| (|1978h . specifically the B3 curve for 
the B z as tabulated by BC03. We note that the trans- 
formations for the ACS passbands Veo6j 2775 and £850 all 
match closely the transformations from B06. We explic- 
itly note here that we do not attempt to use the BC03 
templates to model the photometry directly. Rather, we 
use the templates to compute a grid of observed and 
redshifted magnitudes to which we fit a simple linear or 
quadratic relation. 

The t models used for the transformations span a range 
of parameters. These models had exponential time-scales 
of 0.1 to 5 Gyr, ages from 0.5 Gyr to 12 Gyr and three 
metal abundances, 2.5 solar, solar and 0.4 solar. We 
then fit a linear or quadratic relation between the ob- 
served color or magnitude and the redshifted magnitude 
or color, for all models including all ages of those models. 
For the ACS transformations, we restricted the transfor- 
mations to the range in colors we observed for galaxies 
at the redshift of the cluster and tabulate them in Table 
[TJ All observed magnitudes use the AB system unless 
otherwise noted. 

2.3. Spectral Observations 

The spectral data were acquired over two observing 
runs in 2008 January and March, in excellent seeing and 
generally photometric conditions. The exposures were 
20 minutes long using the 600 line mm -1 grating on the 
DEIMOS spectrograph. The slits were 1" wide, yielding 
a resolution of ^3.8 A or an instrumental a =60 km s _1 . 
The slits were at least 8" long, and between each expo- 
sure we offset the telescope by 2" along the slit. There- 
fore, the object does not land on the same detector pixels 
in sequential exposures, but is instead offset by 2". Each 
mask was observed for six exposures, or w hr. The final 
exposure was 16 hr for the deepest galaxies, 10 hr for the 
brighter objects. 

2.3.1. Spectral Data Reductions 

Our sky subtraction and spectral extraction used the 
techniques outlined in Kelson (2003) as implemented by 
D. D. Kelson. We used the internal quartz lamps to flat 
field the data. Wavelength solutions for each observa- 
tion were determined using the night sk y lines. Each 
slitlet has an independent sky model, see [Kelson] ([2003D 
for details, which is then subtracted from each 20 minute 
observation. 

Bright stars in each mask are used to measure the spa- 
tial offset between frames, as we offset the telescope be- 
tween each exposure. After the sky model was subtracted 
from each frame, we then subtracted sequential frames to 
remove any residual sky signal. Each spectrum in each 
frame has a separate model used for optimal extraction 
(|Hornelll986l) . The initial guesses for the profiles of the 
spectra are based on the bright stars in each mask, but 
are then fit to each spectrum in each exposure. The final 
one dimensional spectrum is an extraction based on all of 
the separate exposures simultaneously, re-binned to the 
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Table 1 

Transformations from Observed to Redshifted Magnitudes 



Filter a 


Transformation 


Color Range 


(mag) 


(mag AB) 


(mag AB) 


r z 


= J + 0.118(7 — J) + 1.538 b 




B z 


= i 775 - 0.113(V 6 06 " 2 8 50) + 0.827 


(^606 - 2850 ) < 1-5 


B z 


= i 775 - 0.190(V 6 06 - zsso) + 0.964 


1-5 < (V 6 06 - 2 8 50) < 2.6 


(U-V) z 


= 0.876(V 6 06 - 2850) - 0.757 


(Vk>6 - Z850) < 1-5 


(U-V)z 


= 1.002(^606 " 2850) - 1-063 


1-5 < (V 6 06 -2850) < 2.6 


B z 


= / 8 i4 - 0.003(V 6 06 - /si4) + 0.778 


V606 — -^814 < 1-2 


B z 


= / 8 i4 - 0.006(y 6 06 - ^814) + 0.783 


1.2 < V 60 6 - hu < 2.4 


(U-V) z 


= 1.135(V 6 06 - hn) - 0.861 


V(i()(i — Isi4 < 1.2 


(U — V) z 


= 1.428(V 6 o6 - ^814) - 1-363 


1.2 < V 60 6 - ha < 2.4 


a Results are in 


Vega magnitudes for the UBV transformations, and we used the Buscr ( 1973) filter curves. 


For the r z , the magnitudes 


are AB and use 


the SDSS filter curves. 


b The observed 


magnitudes are in the Vega system, unlike the other photometry in this table. 





typical dispersion of the data, 0.64 A pixel 1 . We plot 
example spectra in Appendix A, Figure Q] 

2.3.2. Measuring the Velocity Dispersions 

Velocity dispersions and high-precision redshifts are 
measured with a direct fitting method (see, e.g., 
iKelson et al.l l2000bh . as described by Ivan der Wei et al.l 
( 2004f ) and vdW05. using a high-resolution solar spec- 
trum that is smoothed and re-binned to match the galaxy 
spectra. If the continuum of the template and galaxy 
spectra are well matched, template mismatches from dif- 
ferent spectral types introduce errors of, 2% -5%. This, 
combined with a small systematic uncertainty in the 
wavelength calibration and the possible non-Gaussian 
velocity distribution of the stars, prompts us to add 
in quadrature an empirically determined 5% error (see 
vdW05 for a fuller explanation) to the formal fitting er- 
ror on the measured velocity dispersion. Fitting in real 
space, as opposed to Fourier space, has the advantage 
that pixels can be weighted by the signal-to-noise (S/N) 
ratio. In addition, we always mask regions that are af- 
fected by the atmospheric A and B absorption bands 
and emission lines. By default, we also mask B aimer 
absorption features, as these are unsuitable for velocity 
dispersion measurements. However, the He line, which 
is blended with the Ca H line, is not masked if Balmer 
absorption lines are weak, as is the case for the majority 
of the sample for which we managed to measure velocity 
dispersions (26 out of 39, including the three from Wuyts 
et al. 2004 that were re-observed). For the remaining 13 
galaxies, Balmer lines are the dominant absorption fea- 
tures, such that we are forced to mask the Ca H and 
He blended line as well. The final velocity dispersions 
include a 6.6% aperture correction to give the central ve- 
locity dispersion within a fixed physical aperture of 1.6 
kpc or 3.4" at the redshift of Coma, the same aperture 
used by JFK96. 

We note that thanks to the excellent quality of the 
spectra, only two galaxies out of the 55 targeted cluster 
galaxies with ACS imaging do not have measured veloc- 
ity dispersions because of too low S/N in the continuum. 
For a further 14 galaxies we cannot determine the veloc- 
ity dispersion due to the lack of metal features, as their 
continuua are either featureless (with emission lines) or 
dominated by strong Balmer absorption lines. These are 



all classified as late- type galaxies, in agreement with their 
spectral characteristics. Finally, one early-type galaxy 
has no measured velocity dispersion because of a mod- 
erate active galactic nucleus contribution. Because of 
this high success rate, we now have, for the first time, a 
high-redshift sample of early-type galaxies which spans a 
factor of ~4 in velocity dispersion and 2 orders of mag- 
nitude in dynamical mass. All galaxies with spectra are 
listed in Tabled those without measured dispersions are 
so marked. 

2.4. Low-redshift Comparison Sample 

Our low redshift c omparison sample is a subset of the 
Coma sample from iHolden et all ([2007). Each galaxy 
has surface brightness profiles fit to the Sloan Digital 
Sky Survey (SDSS) imaging data in the same manner 
as we have for MS 1054-03. Each galaxy is a member 
of the JFK96 E and SO sample used in that FP study. 
We use the dispe r sions from JFK96 and the colors from 
lEisenhardt et al.l ([2007). The colors from Coma match 
the iBower et all |T992) colors well, but cover a larger 
number of galaxies. We note here tha t we use only the 
half-li ght radius colors from Table 8 of lEisenhard t et al.l 
(|2007fh For the total magnitudes, we use our values from 
the surface brightness profile fits to the SDSS images. 
The larger sample of sizes and colors allows us to expand 
on the sample of JFK96 which had B data for only a sub- 
set of all of the Coma galaxies with dispersions. We list 
in Table @] the magnitudes, rest-frame colors and disper- 
sions we use. Not all of th e samp le of JFK96 have colors 
listed in lEisenhardt et al.l ([2007) as noted in Table |U 

We use the same procedure and t model population 
templates from BC03 as we used for the MS 1054-03, 
see Section 12.2.21 to transform the z = 0.023 observed 
U and V into z — values U z and V z . We find that 
there is some curvature to the relation, and we require a 
quadratic fit of 

(U-V) z = (U -V)-0.067(C7 -V) 2 +0.220(Z7 -F)-0.220 

to transform the observed magnitudes into the redshifted 
ones. For galaxies on the red-sequence of Coma, we find 
U-V ) z = (U-V) -0.04 mag. Using the lColeman et all 
19801) elliptical model, we find (U-V) z = (U-V)-OM 
mag, in good agreement with our quadratic relation. We 
note that these values are significantly different from the 
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{U-V) z = (U-V)-O.OO as found bv fBower et all (fl992T) . 
For the rest of the analysis, we will use the second order 
equation above relating the observed U — V and the z — 
[U - V) z color. 

3. DATA ANALYSIS 
3.1. The Colors of the Galaxy Population 

The colors of the galaxy population in MS 1054-03 have 
been studied extensively in B06. We show in Figure [5] 
the color-magnitude diagram of the population of galax- 
ies that we targeted for dispersions in this study and the 
larger sample of galaxies that could have been targeted, 
a superset of the BOG sample. We note here that we 
remeasured the colors for all galaxies in our sample, in- 
stead of using the values from B06 and we will use these 
colors for our later measurements of the evolution of the 
cluster population. 

Because two different programs observed MS 1054-03 
which used two different sets of filters, we convert the 
colors and magnitudes of galaxies observed in 7gi4 to 
2775 using the following relations 

i 7 75 = ha + 0.248(V 60 6 - hu) ~ 0.211 

V^oe - «775 - 0.752(y 606 - hli) + 0.211 

, which are only valid between 1.2 < Vqoq — I^u < 2.4. 
In Figure [21 we also removed galaxies from our sample 
that do not have spectroscopic redshifts. 

Our colors significant a strong color-magnitude relation 
with a small scatter. The scatter for the whole sample 
of E and SO galaxies with redshifts is a = 0.077 ± 0.008, 
larger than the scatter measured in B06. This larger scat- 
ter comes from the subsample of galaxies with I4o6 — ^814 
colors that have been transformed into Veo6 — *775 colors, 
removing them brings the scatter down to the size as 
measured in B06. This is consistent with the ~0. 01-0. 02 
mag error typically found in these color transformations. 
We find we recover the slope and intercept from B06 with 
our larger sample, showing consistency between the two 
sets of measurements. 

3.2. Success Rate 

Evaluating the degree of evolution in the galaxy popu- 
lation requires calculating the completeness of our sam- 
ple. Part of this is estimating how successful we are mea- 
suring dispersions as a function of galaxy property. We 
targeted a total of 62 galaxies, with 44 above of magni- 
tude limit of J < 21.2. 36 have measured dispersions, 
with 31 above J < 21.2. 

The simplest way of examining the completeness is to 
plot the distribution of galaxies with measured disper- 
sions as a function size, color, and magnitude. First we 
show the color-magnitude plot in Figure followed by 
the size-magnitude relation in Figure [4] In both cases, 
the magnitude is the J, the selection passband. 

3.2.1. The Success Rate of Measuring Dispersions for E and 

SO Galaxies 

We compute the fraction of galaxies with measured 
dispersions as a function of the magnitude and offset 
from the color-magnitude relation, shown in Figure [3l 
Our success in measuring velocity dispersions for E and 



SO galaxies was high, with 29 dispersions out of the 34 
galaxies targeted. The rate of success for measuring dis- 
persions ranges from 100% at the brightest magnitudes to 
70% for the galaxies below our magnitude limit J = 21.2, 
averaging ^80% above our magnitude limit. When we 
look at galaxies as a function of I instead of J magnitude, 
we find a similar distribution. The fraction of galaxies 
targeted with measured dispersions is flat for I < 22 at 
93%, dropping to 70% for the faintest 22.5 < I < 23. 

For E and SO galaxies, we find no evidence of a size 
dependence on the success rate of measuring dispersions. 
This likely implies that magnitude is the important vari- 
able in decidin g the success of a di s persio n measurement. 
Both T05 and Ivan der Wei et all ((20081 ) found that the 
signal to noise of a spectrum was predicted by the scaling 
of FP variables (I e )rl- 9 , which is very close to magnitude 
of a galaxy, {I e )r 2 e . 

3.2.2. The Success Rate of Dispersions for Early-type Spiral 

Galaxies 

Our success with late-type galaxies is lower, with 10 
measured dispersions out of the 20 targeted galaxies with 
J < 21.2, or 50% on average. We find that the success 
rate is flat with size and Sersic index. It is also flat with 
magnitude, but only for galaxies J < 21.2. We measure 
no dispersions for the spiral galaxies below our J < 21.2 
magnitude limit. The most important variable is the 
color, as we are biased toward measuring the dispersions 
of red galaxies. Early-type spirals on the red-sequence 
have a 100% success rate (see Figure [2] for our defini- 
tion of the red-sequence), while blue- ward of that, the 
rate drops to ~30%. Because of this strong color depen- 
dence, we will give bluer spirals galaxies higher weight in 
our later analysis. Measuring the velocity dispersion re- 
quires the presence of strong metal lines, so we are biased 
against young or low metallicity clusters in our sample. 
The small spi ral fraction in MS 1054-03 for moderate 
mass galaxies ([Holden et al .1120071 ) means the overall im- 
pact of this bias will be small. 

3.3. How Representative is the Sample? 

Thirty percent of the galaxies with J < 21.2 that could 
have been targeted had slits placed on them and 26% 
of galaxies with J < 21.2 have measured dispersions. 
These numbers, of course, depend mildly on the apparent 
magnitude of the galaxy. 

The second sample we would like to consider is the 
published work of W04. W04 had a selection limit of 
1 — 22 but preferentially targeted brighter galaxies, with 
more than half of the sample above I < 21 (or J < 19.8). 
This results in a sample with a much brighter effective 
magnitude limit. As can be seen in Figure 02 most of 
the galaxies that could have been targeted, and the ones 
that had been targeted, are on the red-sequence. 

Considering only E and SO galaxies, our spectroscopic 
sample, combined with that of W04, targeted 38% of the 
galaxy population with J < 21.2, with 36% of E and 
SO galaxies above J = 21.2 having dispersions. This 
combined sample of velocity dispersions represents a sig- 
nificant fraction of the E and SO population in MS 1054- 
03. For the rest of the paper, we will use the complete- 
ness fraction of galaxies as a function of J magnitude as 
weights. The inclusion of the W04 sample makes these 
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Figure 3. Color-magnitude relation for the parent sample of galaxies in MS 1054-03. We plot the I — J color of galaxies as a function 
J magnitude. The solid diamonds are the galaxies with measured dispersions from this paper while open circles with "+" symbols are 
galaxies we targeted but did not measure a dispersion for. Solid circles are galaxies with dispersions from W04. We show, as a vertical 
line, our magnitude limit of J = 21.2. The other line is the color-magnitude relation fit to the E/SO population. Along the top, we show 
histograms of the fraction of galaxies with dispersions (filled gray) and targeted (open) as a function of J magnitude. This includes the 
galaxies from W04 - or the solid circles - so this plot shows the relative weight of the different galaxies in our analysis. Along the right side, 
we plot the fraction of galaxies with dispersions and targeted with J < 21.2 by color. We find that our sample of galaxies with dispersions 
covers most of the range of colors and magnitudes in the parent sample. We successfully measure dispersions for ~80% of E/SO galaxies 
and ~50% of spiral galaxies above our J = 21.2 magnitude limit that we targeted. Thus, 26% of all cluster galaxies in the ACS imaging 
and with J < 21.2 have dispersions, matching the visual impression of this figure. For 14 galaxies we targeted, we cannot determine the 
velocity dispersion due to the lack of metal features, as their spectra are either featureless (with emission lines) or dominated by strong 
Balmer absorption lines, producing a decreasing fraction of galaxies with bluer I — J colors. 



weights a strong function of magnitude, as can be seen in 
Figured] These completeness fractions range from 100% 
at the brightest magnitudes to 20% at our completeness 
limit of J = 21.2. 

4. EVOLUTION OF THE FUNDAMENTAL PLANE AND 
GALAXY PROPERTIES 

4.1. Fundamental Plane Evolution 

4.1.1. Measuring the Tilt of the Fundamental Plane 

We compute the parameters of the FP by fitting the 
J < 21.2 magnitude limited subset of our dispersions 
(a), effective radii (r e ) and our B z surface brightnesses. 
To fit the FP, log r e = a log a + /3 log(7 e ) + 7, where (I e ) 
is — 0.4/^b, we minimized the mean absolute orthogonal 
deviation, as was done in JFK96: 

log v e ol log o~ j3 log( I e ) - 7 



A = 



/3 2 ) 



around the FP in all three projections, r e , <J and (I e ) . We 
computed the average values of each of the coefficients 



for our final estimates of a, ft and 7. We estimated the 
errors on all of our fits to the FP by bootstrapping the 
data. 

To compute the best fitting plane, we need include 
our incompleteness. We did this two different ways. 
First, we use our weights we computed above in Sec- 
tion 13.31 Thus, galaxies at the fainter end of the lu- 
minosity function, closer to our limiting magnitude, get 
a higher weight. Second, we trim the data at a fixed 
value of a. Because galaxies become intrinsically fainter 
at lower masses, by trimming at a fixed cr, we remove 
the very lowest M / L galaxies where our sample will not 
have the corresponding high M/L galaxies at the same 
dispersion. We find that at J = 21.2, the typical dis- 
persion is log 10 a = 2.10 ± 0.10, and thus we select limit 
of log 10 a — 2.20. We will examine how our best fitting 
plane varies depending on whether or not we remove the 
lowest a galaxies or use our weights. 

Our resulting best fit values are summarized in Table 
[21 and, for comparison, we give the relations from JFK96 
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Figure 4. Size-magnitude relation for the parent sample of galaxies in MS 1054-03. We plot the J magnitude of galaxies as a function 
HST/ ACS measured half-light radii using the same symbols as Figure [3] The fraction of galaxies with dispersions and J < 21.2 is plotted 
as a function of size in the gray histogram above the plot, and the fraction of galaxies targeted for dispersion measurements is shown with 
the open histogram. We find that fraction of galaxies with dispersions is relatively flat with size above our magnitude limit, with exception 
of the two brightest galaxies in the cluster and the smallest galaxies which are predominately below are magnitude limit. This means that 
we have no strong surface brightness selection bias. This is in accordance with previous work which finds that the S/N of a spectrum is 
proportional to (I e )r\-®, which is close to (I e )r^ or the observed magnitude (T05;[van dcr Wcl ct al .1 120081) . 



for the B selected sample. We verified that we recover 
the FP values from JFK96 when we fit to the values listed 
in Table H 

When fitting the FP in MS 1054-03, at first we re- 
stricted the sample to only those galaxies that are E and 
SO galaxies, a similar selection as JFK96, but using the 
log 10 a > 2.2 and using the weights from above, we find 
a = 1.25 ± 0.11 and ft = -0.76 ± 0.03. The value of a is 
in good agreement with that of JFK96 but the value of ft 
appears different a t the level of two standard deviations. 
iHyde fc Bernard! (|2009f ) show that sample selection can 
cause large changes in the best fitting slopes of the FP, 
and specifically trimming at a fixed a. The JFK96 sam- 
ple is complete to a fainter r limit than the equivalent 
J limit for our sample and did not have a limiting a as 
we did. So, we restricted the sample of Coma galaxies to 
those galaxies that would match our high redshift selec- 
tion criteria, namely E and SO galaxies with an equivalent 
r selection corresponding to our J < 21.2 at z = 0.83 and 
the same limiting a. We refit the FP for the Coma sam- 
ple, and we find that value of ft for the Coma galaxies 
decreases to (3 — —0.78 ± 0.04 with a unchanged. Thus, 
we conclude that the mild difference in ft we find in com- 
parison to JFK96 is likely a result of sample selection, 



and that there is no evolution in the tilt of the FP. 

Including early-type spirals in our sample does not 
change the result signifi c antly, a s was seen in a z = 0.33 
cluster bv IKelson et all (|2000d) . Fitting the FP to all 
galaxies with log 10 er > 2.2 and using the weights from 
above, we find a = 1.18 ± 0.16 and ft = -0.76 ± 0.06. 
Lowering the limiting a to log 10 a > 2.1 and using 
the weights from above, we find a — 1.19 ± 0.13 and 
ft = — 0.78 ± 0.05. In general, these numbers are in good 
agreement with the values for a more restrictive subsam- 
ple of E and SO galaxies with log 10 a > 2.2. This shows 
that our results are robust to different limiting velocity 
dispersions and uncertainties in the sample selections. 

From these tests we conclude that the observations are 
consistent with an unchanging tilt in the FP. In Figure 
[5l we plot the FP for our z = 0.83 sample. We use the 
values of a and ft from JFK96 to plot both the low and 
high redshift samples. 

4.1.2. Residuals Around the Fundamental Plane 

vv07 found that there was a strong correlation be- 
tween the residuals around log r e and other quantities 
related to the mass of the galaxy, such as a or the dy- 
namical mass. Such a deviation could be caused by cur- 
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Table 2 

Summary of FP Parameters 



Cluster Limiting logo - Morphology a ft 

(km s _1 ) 



JFK96 a E and SO 1.20 ± 0.07 -0.83 ± 0.02 

Coma b 2_2 E and SO 1.18 ±0.08 -0.78 ± 0.04 

MS 1054-03 2~2 ATI 1.18 ±0.16 -0.76 ± 0.06 

MS 1054-03 2.2 E and SO 1.25 ±0.11 -0.76 ±0.03 

MS 1054-03 2.1 All 1.19 ±0.13 -0.78 ± 0.05 



a The best fitting slopes for the B from the r selected Coma sample from Table 5 of JFK96. 

b The best fitting slopes as determined by for the Coma sample from Table [4] after selecting galaxies above the log a limit and above the 
equivalent r magnitude limit that matches our J magnitude limit at z = 0.83. 



vature in the FP and is more thoro ughly explored in 
Ivan der Marel fc van Dokkuml (2007a b). We assume the 
FP relation from JFK96 and we plot, in Figure [51 the 
deviation between the measured r e and the predicted r e 
from the FP as a function of r e , er, and the dynamical 
mass Mdyn- We computed Md yn is in M using the re- 
lation Mdyn = 5a 2 r e /G or 

log 10 M = 2 log 10 a + log 1Q r e + 6.07 

with a is in km s _1 and r e is kpc. We find galaxies with 
dynamical masses Md yn < 10 11 Af Q , that lie close to the 
z = relation after adjusting for passive M/L evolution 
(see below). This result is expected given the lack of 
evolution we find in the tilt of the FP. 

4.1.3. M/L Evolution from the Offset of the Fundamental 

Plane 

Traditionally, the offset in the FP is used to measure 
the amount of M/L evolution. This is done by assuming 
that the offset, the value of Alog 10 r e , is all caused by 
luminosity evolution, or a change in (I e )- If we take the 
value we measure, Alog 10 r e — — 0.38 ± 0.02, in Figure|5] 
we find Alog 10 M/L B = -0.44±0.03 or d\og 10 M/L B = 
—0.60 ± 0.04 dz. Our values are in excellent agreement 
with the value of A log 10 Af/is = -0.415 ± 0.040 at 
z = 0.83 as found by vv07 but much smaller than the 
dlog 10 M/L B = -0.72toi 7 5 dz of T05. 

4.2. Evolution in the M/L from the Virial Estimator 

We can combine our measurements of the total B light 
with the M to estimate a M/ Lb value or 

M/L B = 5a 2 /(2nGr e (L e )) 

We show in Figure [7] the relation between M /L b and 
a for Coma and our sample of galaxies in MS 1054- 
03 from the virial M/L estimator. We fit the re- 
sults at both low and high redshifts using a robust lin- 
ear fitting technique that minimizes the median abso- 
lute deviation and estimate our errors with bootstrap- 
ping. The bootstrapping is performed in a weighted 
manner, based on the weights determined in Section 
13.31 It appears that the slope has not evolved between 
the low and high redshift samples. In our sample of 
galaxies in MS 1054-03, we find for the Coma sam- 
ple, M/Lb oc (j - 98 ^ - 10 i n excellent agreement with 
the low redshift result of Ivan der Marel fc van Dokkuml 
(|2007bft . while for our sample in MS 1054-03 we find 



M/Lb oc <jL12±o.21_ if we assume no change in the 
slope, we find Alog 10 M/Lb = —0.50 ±0.02 at a fixed a 
for the whole of the galaxy population with J < 21.2. 
For the red-sequence, E and SO population, we find 
Alog 10 Af/iB = — 0.50 ± 0.03, or no measurable differ- 
ence. 

Our measured M/Lb evolution using the virial estima- 
tor is higher than the Alog 10 M/L B = -0.44 ± 0.03 we 
find in Section 14.11 In principle, these results should 
be the same if the evol ution in the FP is caused en- 
tirely b y M/L evolution. Ivan der Marel fc van Dokkuml 
(2007b) showed that structural evolution, such as caused 
by higher rotation rates or size evolution, would cause 
the disagreement we obs erve between the two estima- 
tors of M/Lb evolution. Ivan der Marel fc van Dokkuml 
(2007b) find, for their sample, that the virial estimator 
yields a result closer to t he M/L as estimate by the dy- 
namica l modeling done bv lvan der Marel fc van Dokkuml 
(2007tJ). For the rest of this paper, we will use the virial 
estimator as a consequence, but will note how using the 
FP results changes our results. This s hould avoid some 
of t he bias found inlSaglia et al.l (120101). Finally, we add 
that Ivan der Wei fc van der Marell (|2008D found little ev- 
idence for higher rotation rates in field E and SO galaxies, 
thus evolution in rotation alone does not explain the re- 
sults of T05 or vdW05. 

4.3. Color Evolution 

In Figure [U we show that, as expected, the rest- frame 
(U - V) z colors of the MS 1054-03 galaxies are bluer 
than the colors of the Coma galaxies. This shift can be 
seen both at fixed luminosity and at fixed a. Since a 
is expected to evolve little or not at all as compared to 
luminosity, we adopt the color evolution at fixed a as the 
quantify of interest. For Coma, we find the relation (U — 
V) z = 0.653±0.093 log 10 (cr/200 km S - 1 ) + 1.475±0.015. 

We find no ev olution in the slo pe, as expected from the 
larger survey of iMei et al.l (|2009[) which found no change 
in the slope of the color-magnitude relation for eight z > 
0.8 clusters. We measure the offset in the average color 
for all red-sequence E and SO galaxies in our sample with 
log 10 (cr) > 2.2 and we find A(U - V) z = -0.24 ± 0.02 
mag between z = 0.831 and z = 0.023. As above, we 
use bootstrapping to estimate the errors on the color-cr 
relation in Coma and on the amount of color evolution. 

5. THE EPOCH OF CLUSTER GALAXY FORMATION 

We use a combination of the M/Lb evolution and the 
U — V color evolution to estimate the median epoch of 
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Figure 5. Fundamental plane for the Coma sample of JFK96 (blue open squares) and our sample in MS 1054-03. Spirals are marked by 
green spiral symbols while E and SO galaxies are filled red points, circles for those on the red-sequence and triangles for those off. We plot 
all the galaxies in Coma, regardless of a or r magnitude. We use the values of the slope of the FP, a = 1.20 and = 0.83, from JFK96, 
and plot the FP relation as a blue solid line for Coma. The red solid line is the FP relation for MS 1054-03 using the same slope as we 
do for Coma, but with an offset to match the relation in MS 1054-03. We find no evidence that the tilt in the FP evolves, as we recover 
the relation of JFK96 for the galaxies in our sample. We do find that the FP has a different zero-point. The change in the zero-point 
is Alog 10 r e = —0.38 ± 0.02, which corresponds to evolution of Alog 10 M/Lg = —0.44 ± 0.03 for the stellar populations of the cluster 
galaxies. 



star-formation for galaxies in our cluster sample. We es- 
timate the luminosity weighted age of the stellar popula- 
tions by c o mpu ting which single stellar population from 
iMarastonl (|2005l hereafter M05) best reproduces both 
M/Lb and U — V evolution we observe. We will assume 
that the populations are coeval and we include all galax- 
ies, regardless of morphological type, in MS 1054-03. Fit- 
ting both the M/Lb evolution and the color-evolution 
simultaneously, we find that M05 models with a solar 
abundance prefer z* = 1.81q; 2 . Because of our complete- 
ness limits, this is in effect the median age of formation 
for galaxies with log 10 a > 2.2 or, roughly, a dynamical 



mass of M > 6 x 10 10 M Q . 

If we use the smaller Mj L from the offset in the FP, 
see Section 14.11 we find a higher formation redshift, 
= 2.0lg'2 > though well within the range of the er- 
rors. In fact, the uncertainties in this are as much from 
the assumptions about the stellar populations as from 
the errors on the measurements. For example, assum- 
ing a metallicity logZ/Z© = 0.35 increases the forma- 
tion epoch to = 2.5^Q 2 • We note that the pace 
of M/L evolution and the pace of the U — V evolution 
are in mild disagreement with expectations of a higher 
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Figure 6. Deviations from the FP of JFK96 as a function of r e , a, and mass (M^ yn = 5a 2 r e /G). We use the same symbols as in 
Figure \5\ and show the approximate J = 21.2 selection limits for a red-sequence galaxy at z = 0.83 with dashed lines. Galaxies below 



those lines are in our fainter subsample (21.2 < J < 21.8). The average offset of the FP, Alogj^ 



-0.38 ± 0.02 (which corresponds 



to A log 10 M/Lb = —0.44 ± 0.03), is illustrated with a red line. The offset of z = 0.83 galaxies from the z = 0.023 FP is not larger for 
galaxies with low a or mass as long as they lie above our limiting magnitude, as expected if the slope or tilt of the FP does not evolve. 



metallicity population, but if we assume the slower pace 
of evolution from the offset of the FP, we find a better 
agreement, though we prefer the estimate from the virial 
M/L estimator for reasons discussed in Section |4~2"1 

There are significant differences between the models 
of BC03 and those of M05, but these are in how the 
models address the post-main-sequence evolution, specif- 
ically the thermally pulsating asymptotic giant branch 
stars. This can cause significa nt differences in the near 
IR fluxes, but as was shown in Ivan der Wei et "ail (2006) 
the impact on the M05 predictions for the optical prop- 
erties of older stellar populations is small. We estimated 
the typical epoch of formation for the whole sample us- 
ing the BC03 models and found z* = 1.6^"'^ , in good 

agreement with the value of z* = 1.8^ ' 2 found using 
M05 with the same assumptions. 

From the above results we conclude that the epoch 
of formation for the galaxies in MS 1054-03 is z* — 
2.0 ± 0.3 ± 0.3 (sys), in good agreement with the results 
from vv07 for their sample of higher mass (> 10 11 M©) 
galaxies. The systematic error shows the range of al- 
lowed values given both the uncertainties in the stellar 
population models as well as the possible different M/L 
values. 

5.1. Constraints on the Initial Mass Function 

The evolution of M/L for passively evolving systems 
depends critically on the IM F, while the ev olution in op- 
tical colors depends less so (Tinslcv Il972h . This comes 
about because both the colors and luminosity of a pas- 
sively evolving stellar population are determined by the 
properties of the stars at the main-sequence turn-off. The 
color evolution depends mostly on the effective tempera- 
ture of those stars, and thus is most strongly influenced 
by the age and typical metallicity. The luminosity evo- 
lution, however, depends on not just luminosity of the 
individual stars, but the relative numbers of stars as a 
function of mass. Thus, our measurements of the ratio 
of the pace of color and luminosity evolution can be used 
to constrain the IMF of stars at the main-sequence turn- 
off, around 0.8 — 1.0M Q for old stellar populations in our 



target galaxies. 

vD08 found that, for cluster galaxies with masses > 
10 11 M Q , the color and M/L evolution favored an IMF 
very different from a Salpeter slope in the mass range of 
0.8 — l.QM©. The assumptions we made in the previous 
section use a Salpeter-like IMF to estimate the formation 
epoch from the M/L evolution. vD08, however, shows 
that the IMF that most closely agrees with the sample 
in that paper yields a formation epoch of z* = 3.7^q g 

for cluster galaxies with > 10 11 M©. In light of the 
results of vD08, we have examined our data removing 
the assumption of a fixed IMF. 

vD08 provides a useful relation for constraining the 
IMF with the evolution in the color and M/Lb- vD08 
finds that for solar metal abundances with the M05 mod- 
els with ages >1 Gyr 

, where the terms Alog 10 (M/Ls) and A(U — V) rep- 
resent the amount of evolution in the M/Lb and color, 
res pectively. The t erm x is the slope of the IMF, where 
the Salpeter (|1955D value is x = 1.35. 

We plot Alog 10 M/Lb and A(U - V) along with the 
expected evolution in Alog 10 M/ Lb and A(?7 — V) in 
Figure [9] We find that our data are in mild disagree- 
ment with the expected evolution from a Salpeter IMF, 
with a slope of IMF better described by x = 0.9 ± 0.2, 
shown in Figure |9] with a red line. This is a difference 
of 2.25 standard deviations from the Salpeter value of 
x = —1.35. This value modestly changes the best fitting 
formation epoch from = 1.8lg 2 *o z* = 2.0 ±0.3. As 
we discuss above, such a shift is in line with the other 
systematic uncertainties. 

5.1.1. Contrasting Our Results with vD08 

vD08 found much lower values for x, with a preferred 
value of x = —0.1 for galaxies with a solar metallicity. 
We find x = 0.9 ± 0.2 in large part because we compare 
only Coma and MS 1054-03. If we include the lower red- 
shift data from vD08, specifically the data from A2218, 
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Figure 7. M/Lb as measured by the virial estimator (M/Lg = 5a 2 / (2nGr e (I e ))) vs. a for the samples of Coma and MS 1054-03. We 
use the same symbols as in Figure [5] The black dashed line is the approximate selection limit for a red-sequence galaxy in MS 1054-03. 
We show the typical errors for MS 1054-03 in the lower right. Our high redshift slope is compatible, within the errors, with the slope we 
find at low redshift. We find a relation of log 10 M/Lg = 0.981og 10 (a-/200 km s —1 ) + 0.90 for Coma (shown as a blue line) and offset of 
Alog 10 M/L B = -0.50 ± 0.03 for red E and SO galaxies in MS 1054-03 (shown as a red line) and Alog 10 M/L B = -0.50 ± 0.02 for the 
whole sample of J < 21.2 galaxies, regardless of color or morphological type. 



we would find a shallower slope of x = 0.5±0.2 illustrated 
with a green dash-dotted line in Figure [9] The key differ- 
ence between our result and those of vD08 is not in our 
high redshift data, but rather the lower redshift sample. 

vD08's galaxy samples were selected by including all 
galaxies with dynamical masses > 10 11 Mq. This is in 
contrast to our sample in MS 1054-03, which is limited to 
those galaxies at log 10 a > 2.20, roughly 6 x 10 10 Mq. If 
we instead use the criteria of vD08, namely a subsample 
of E and SO galaxies in MS 1054-03 with > 10 11 Mq, we 
find a AU — V = 0.22 ± 0.04, similar to the AU — V = 
0.20 ± 0.04 found in vD08 and 0.03 mags smaller than 
the offset we measure for the full sample of galaxies in 
MS 1054-03. So part of the difference in our derived slope 



of the IMF comes from our sample probing a larger mass 
range. This sample selection results in a difference in 
AU — V = 0.03 mags or a slope of x = 0.5 ± 0.2 instead 
of x = 0.9 ± 0.2. The rest of the difference between 
our value of x = 0.9 ± 0.2 from our data alone and the 
x = — 0.1 as measured by vD08 is the larger sample of low 
redshift clusters that vD08 used. If we use the same mass 
selection as vD08, all E and SO galaxies with > 10 11 Mq 
and use the low redshift sample of Coma and A2218, we 
recover the x = —0.1 found by vD08. 

The data we have used in this paper have been ana- 
lyzed in a consistent manner throughout. vD08 used a 
larger set of data from t he literature that was made as 
consistent as possible (see Ivan Dokkum fc van der Marell 
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Figure 8. Rcdshifted color-magnitude and color-tr relations. We show the (U — V) z colors both for Coma galaxies (z = 0.0231) and 
galaxies in MS 1054-03 (z = 0.831). We use the same symbols here as in Figure [2] The typical error in the color and a is shown, as a 
function of a, across the bottom of the right panel. We note that the large a sample is dominated by W04 which have larger errors in a 
owing to shorter effective exposure times. We find evolution in the offset of the color-magnitude relation and the color-cr relation, which 
implies a younger population at higher redshift. This evolution is determined by fitting the color-magnitude and color-cr relation to the 
Coma galaxies, with the relations shown as blue lines, and then computing the best-fitting offset for those galaxies on the red-sequence 
in MS 1054-03. The relations for MS 1054-03 are shown as red lines. We find no statistically significant evidence for a change in the 
slope of either relation when we fit for them to red-sequence cluster members. For all red-sequence E and SO galaxies in our sample with 
lo gio(°") > 2 - 2 ' we find A (U ~ V) z = -0.24 ± 0.02 mag of color evolution. 



[20071 for details). In vD08, it can be seen, both by ex- 
amining the equivalent of Figure [9] from that paper or 
from other measurements of the color and M/L evolu- 
tion, that the low redshift data are as important to the 
conclusions as the high redshift data, a result reinforced 
by the change in the slope of the IMF we find when 
we include the eight galaxies of A2218. Because of this, 
the large samples of low redshift FP samples, such as 
iHyde fc Bernar di (2009), should be used to more fully 
explore joint evolution in M/L and color. From this we 
will have a much more solid measurement of the IMF for 
E and SO systems. 

5.2. Implications from the Lack of Evolution in the Tilt 
of the Fundamental Plane 

We find no evidence for any evolution in the tilt of 
the FP, nor do we find evidence for a change in slope 
of the color-cr relation. Therefore, our data would rule 
out the scenario where the median star-formation epoch 
depends strongly on galaxy mass. We quantified this 
by fitting a model to the color and M/L data in both 
Coma and MS 1054-03. This model assumes that the 
galaxies in Coma and MS 1054-03 is coeval and have the 
same metal abundance, and leaves the scatter as a free 
parameter at each redshift. The model parameterizes the 
formation epoch as power-law function of the a, tf = 
i (a/ 160km s -1 )". Fitting this model to those galaxies 
above our magnitude and a limit, we find a = 0.2+q '\ 5 



and recover a t corresponding to z* = 1.5. The low value 
for a means that, in our best-fitting model, the most 
massive galaxies (log 10 cr — 2.5) have a z± = 2.3+J;|. 
Our resulting value for a is on the small side of those 
reported in the literature, which we will discuss in later, 
in Section EH 

There is an important caveat to this analysis. We 
assume that the whole population in both clusters are 
coeval. Because the cluster population grows by accre- 
tion of galaxies from outside of the cluster, however, the 
population of galaxies in MS 1054-03 represents only a 
subset of the z = cluster galaxy population. Thus, 
we would overestimate the epoch of cluster galaxy for- 
mation as we would only include those galaxies that are 
already E, SO or early-type spirals at z = 0.83. If this 
process is dramatic, we would expect to find rapid evolu- 
tion in the fraction E and SO galaxies in clusters. How- 
ever, for mass selected samples, the amount of observed 
evolution in the fraction o f E and SO galaxies is minima l 
(IHolden et al]l2?Ml2"fi07h . Ivan Dokkum fc Franxl (pOOlh 
provide a method for calculating the impact of this bias. 
We find if we assume that 90% of the galaxies were in 
our sample at z = 1, there is no shift in the epoch of 
formation. Assuming that fraction is only 50%, the ef- 
fective epoch of formation shifts downward by Az = 0.2, 
= 1.4 for galaxies with log a — 2.2. It is likely that this 
fraction actually varies with velocity dispersion, so, in ef- 
fect, the actual slope of the relation between formation 
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Figure 9. Alog 10 M/Lg vs. A(C7 — V) z for our Coma E and SO sample (blue squares) and for the red-sequence E and SO galaxies in 
MS 1054-03 with log 10 a > 2.2 (red, open circles). For each galaxy, we measure the offset with respect to the Coma M/Lg — a (Figure 
[TJl and U — V-cf (Figure [SJ relations. Our mean offsets in A(£/ — V) z and AM/Lg with ler error in red and 3 a in lighter pink, with our 
corresponding 1 and 3 cr error ellipses for Coma in dark and light blue respectively. The typical individual errors for MS 1054-03 are shown 
by a red error-bar in the lower right. We show the expected evolution in Alog 10 M/Lg vs. A(U — V) z for an IMF with a Salpeter slope 
at ~ 0.8 — I.OMq with a metallicity ranging from solar to more than twice solar, log 10 (Z / 'Zq) = 0.35, in the filled gray region. The best 
fitting IMF slope, x = 0.9 ± 0.2 from only our data is shown by the red line. Our sample, alone, shows only tentative evidence of a shift in 
the IMF, at the level of 2.25 standard deviations. The black error bars are the data from vD08, with the expected evolution from the best 
fitting IMF of vD08, x = —0.1, shown as a dashed line. Using the lower redshift data set of vD08, namely Coma and A2218, we would find 
x = 0.5 ± 0.2, shown as a green dash-dotted line. The general agreement with our results (in color) and those of vD08 (in black), however, 
implies that both samples are drawn from similar parent samples. From this, we conclude that the stellar populations should have similar 
IMF slopes. 



epoch and dispersion is mildly steeper than we measure. 
If we assume that 50% of the logo- — 2.2 galaxies were 
not present in the sample while 90% of the loger = 2.5 
were present, this shifts the resulting a from a = 0.2^q '\ 5 



to a — 0.25^q^ 5 , a very modest change. 

5.2.1. Implications for the IMF 

iRenzinil (|2006f ) shows that a trend in the IMF with 
mass would cause the tilt in the FP to evolve with red- 
shift. Because our measured slope of the FP changes 
little, if at all, we can immediately conclude that the 
IMF of the lower mass galaxies in our sample is similar 
to that of the higher mass systems. In Section 15.11 we 
show that most of the difference between the slope of the 
IMF from our sample, x = 0.9, and the slope of the IMF 
as determined by vD08, x = —0.1, comes from the larger 



sample of clusters used in vD08. The remaining differ- 
ence, a change in the slope Ax = 0.4, comes about from 
our sample of lower mass galaxies in MS 1054-03. Thus, 
the IMF for cluster galaxies has, at most, a mild mass 
dependence for galaxies with log 10 a > 2.2, and a mass 
de pendence that ap pears in agreement with the results 
of lTreu et all (f2010h . 

As many lower mass cluster galaxies must have formed 
in more fiel d-like environments before falling into the 
cluster (eg. iPatel etaD l2009bt iBerrier et alJ I2009L and 
references therein), then the IMF for many field E and 
SO galaxies should be similar to what is found for clus- 
ter galaxies at similar masses. This follows regardless of 
whether the conclusions of vD08 are corre ct or, in con- 
trast, the conclusions of iTreu et al.1 (|2010| ) are correct. 
Therefore, if there is a different IMF for E and SO galax- 
ies, it should be prevalent in all moderate to high mass 
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galaxies regardless of environment. 

5.3. Comparison with Color Measurements 

A number of color-based measurements of the forma- 
tion epoch of cluster galaxies have been made. B06 found 
that the scatter in the colors in MS 1054-03 and another 
z = 0.83 cluster, RX J0152-13, implied a star formation 
epoch of Zj, ~ 2 .2, in good agreement with our results. 
IMei et al~l (|2009Q found a similar scatter in the colors of 
the E and SO population in a sample of eight clusters, 
including the two in B06, indicating a simi lar formation 
epoch for the z ~ 1 cluster population. iHolden et al.l 
(2004) measured z+ = 3+i by fitting the evolution of 
the zero-point of the color-magnitude relation of E and 
SO galaxies, as opposed to the scatter as was done in 
B06, for a sample of 24 clusters spanning a re dshift 
range of 0.3 < z < 1.3. lEisenhardt et al.1 (|2008[ ) and 
iMancone et al.l ()2010[ ) find a similar formation epoch for 
the galaxies on the red- sequence in a much larger sample 
of z > 1 clusters than in lHolden et al.l ([2004D . All of these 
results all point to z = 2 — 3 being the important epoch 
for the star formation of clus t er ear ly- type systems. In- 
terestingly, lEisenhardt et al.l (|2008l ) find evidence that 
the average formation epoch of the cluster red-sequence 
increases for the higher redshift cluster galaxies, the ex- 
pected result if more galaxies are joining the red-sequence 
over time. 

At z ~ 2 — 3, Brammer fc van Dokkuml 

d200l. IZirmet alJ <j200"8h . IKriek et all (120081) 
iBrammer et al.l ( 2003 ) all find a well-defined 



and 
red- 
sequence in both clusters and the field environment. 
These early-type systems have finished forming their 
stars at z ~ 3, a time higher than the mean age as 
measured by colors or by the FP. However, these galaxies 
represent only the oldest 10% of the early-type popula- 
tion, and so are consistent with a more typical epoch of 
z+ ~ 2 for the majority of the population. Therefore, 
color measurements show that ~ L* galaxies have stellar 
populations with formation times of z t ~ 2 — 3, an 
epoch when a sizable fraction of the stars that exist 
today were formed and entirely consistent with our FP 
measurements . 

5.4. Comparison with Archaeological Studies of 
Present-day Galaxies 

ITrager et al.l (2000b. a) found that age and dispersion 
were part of a hyper-plane also including metallicity and 
a enhancement. In general, higher dispersion galaxies 
have older population ages at a fixed metallicity and 
galaxies in denser environments show older ages. A num- 
ber of papers have built on these results, finding a strong 
correlation between the properties of stellar populations 
and velocity dispersions of galaxies for very large sam- 
ples (INelan et al.ll2005HGallazzi et a l]|2006t [Gr aves et al.l 
I2009bllal ; iSmith et al.l I2009D . Interesting ly dynamical 
mass is not th e important variab le, rather [Graves et "all 
(|2009bllah and ISmith et al.1 (|2009l) both found that a is 
the driving parameter for determining the stellar popu- 
lations of z = early-type galaxies. 

The ages derived from these measurements are strik- 
ingly young for lower dispersion galaxies, roughly 30%- 
50% of galaxies at log 10 a = 2.2 having form ation ages 
corresponding toz~ 0.8. ITrager et al.l (J2008) found sim- 
ilar results for E and SO galaxies in the Coma cluster, a 



typical age of 5-8 Gyr, or roughly a z* ~ 0.5 — 1, re- 
gardless of galaxy mass. Correspondingly, many authors 
find a steep relation between the age of the stellar pop- 
ulation and velocity di spersion, with logar i thmic slopes 
of 0. 4 - 0.6 common (ITrager et al.l I2000at INelan et all 
120051: ISmith et al.l I2009 D and values up to ~ 1.1 found 
( Bernardi et al.l l2003bO . Kelson et all ([2006f ) found a 
strong relation between metallicity and a, but, in con- 
trast, that after accounting for the selection limits the 
population of CL 1358+62 showed no age variation with 
a. 

Part of the reason for the difference in our results, 
a shallow age-dispersion relation, and the steeper ones 
found by others comes about from our use of broad band 
photometric data as co mpared to fitting mod els of stellar 
populations to spectra. iTortora et al.l (|2009D found that 
the slope of the age-cr relation is much shallower when 
stellar po pulations are fit to o nly broad band photometry 
Further. ICooper et al.l (|2009() found that galaxies in high 
density environments with the same broad band color 
and magnitude as galaxies in low density environments, 
had both larger ages and larger metallicitie s using mod- 
els fit to spectra from IGallazzi et all (|2006l ). This shows 
that the population parameters derived from spectra are 
not th e same as those deriy e d from color measurements 
alone. ITrager fc Somerviliel (|2009l ) provides a potential 
frame work to interpret these seem ingly contradictory re- 
sults. ITrager fc Somerviliel ((2009) find that the ages that 
result from fitting simple stellar population models to 
spectral features are likely when the last 5%-10% of the 
stars formed. Thus these young apparent ages are not 
when the median star was formed, but rather the epoch 
when the last star formation occurred. 

5.5. The Field versus the Cluster Early-type Population 

Our sample consists entirely of cluster galaxies, so 
part of the reason we find a different result from 
field samples s uch as T05 or yd W0 5 c ould simply be 
environ ment. IBernardi et al.1 (l2003al). IBernardi et al.l 
(120061) . ICooper et al.l (|2009f ). lLa Barbera et al.l (|2010f ) 
and lSaglia et al.l ( |2010l ) all find differences in the FP or 
the stellar populations of galaxies such that galaxies in 
clusters appear older than the same mass field galaxies. 

Interestingly, T05 show, in Figure 14 of their paper, 
that most massive (> 10 11 M Q ) galaxies have a forma- 
tion epoch of z ~ 2 — 3. vv07 perform a separate analysis 
and find a similar result. The main differences between 
our results and T05 are at < 10 11 Mq. We find that 
~ L* cluster galaxies have a typical formation epoch of 
= 1.7 to z* = 2, whereas T05 find ~ 1.2 for similar 
mass galaxies in their sample. In our MS 1054-03 sample, 
we have galaxies that show rapid enough M/L evolution 
to be consistent with a formation epoch of z* ~ 1.2, but 
the typical low a galaxy has a much higher M/L (see 
Figure [7|). In the context of our model presented in the 
previous section where t* oc (cr) a , the results of T05 pre- 
fer an«~ 0.4 — 0.45 as opposed to our a ~ 0.2. There 
are two possible explanations for why our results differ 
from T05. First, the selection limits of T05 z band filter 
galaxy selection, as opposed to our redder J, biased the 
resulting sample towards bluer galaxies than our sam- 
ple,. The second possibility is that the results of T05 are 
different for lower mass field galaxies because lower mass 
field galaxies have a different star-formation history as 
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compar ed wi th lower ma s s clust er galaxies. iPatel et al.l 
(|2009bf) and IPatel et al.l (|2009aD do find that field and 
cluster galaxies have different star-formation rates and 
color distributions at z = 0.83, possibly explaining the 
difference between our results and T05. 

T05 explain the apparently young stellar populations 
seen in lower mass galaxies because of small recent 
episodes of star-formation and not because the galaxy 
po pulation is, as a wh o le, yo ung, akin to the model 
of iTrager fc Somervilld (|2009| ). This picture is given 
m ore support when co mpared with the recent results 
of iThomas et al.l ((2009) , which find that a small frac- 
tion of galaxies in the local universe with an E and SO 
morphology appear to have had recent, small episodes of 
star-for mation. Labeling thes e young galaxies as "rejuve- 
nated" , IThomas et al.l (|2009l ) found the fraction of these 
galaxies increases at lower velocity dispersions, mimick- 
ing the steep a ge-dispersion relatio n found in T05 and 
in other work. iRogers et al.l (|2010[) also found a higher 
fraction of galaxies with recent star-formation in lower 
mass dark matter ha l os. Al l of this is consistent with the 
ITrager &: Somervillel (120091) pict ure . To explain the re - 
sults of. sav. iCooper et al.l (|2009t ) or IRogers et al.l (|2010h . 
these events would be more common in field like environ- 
ments. These are likely the handful of galaxies we ob- 
serve with bluer colors and lower M/Lb values at lower 
values of er. This picture provides a natural explana- 
tion for the build up of the re d-sequence fr o m z ~ 1 
until today, see, for example, lHarker et al.l (|2006| ) or 
iRuhland et al.l (|2009f ) , and could provide th e progenitors 
for th e apparently young Coma galaxies in ITrager et al.l 
(|2008h . 

6. SUMMARY 

We have compiled a sample of velocity dispersion mea- 
surements of cluster galaxies spanning a broad range in 
mass at z — 0.831. Our sample goes much farther down 
the mass function than many previous cluster and field 
samples. Our sample was selected only by J magnitude, 
effectively in rest- frame r at this redshift, closely mimick- 
ing the selection used for many low-redshift samples. No 
morphological or color information were used in the se- 
lection. With HST/ ACS imaging in the bandpass closest 
to the rest-frame B, we determined the size and surface 
brightnesses of the galaxies in our sample. We found that 
our success at measuring velocity dispersions depended 
mostly on the magnitude of the galaxies, with no appar- 
ent bias in size or surface brightness. We did find that, 
for the early-type spirals in our sample, we are less likely 
to measure the dispersions of galaxies with bluer colors, 
so we proceeded by giving those galaxies higher weight 
in our analysis. 

We determined from our sample the best fitting FP 
logr e = alogfj + 61og(/ e ) + c, and measure the color and 
M/L evolution of the early- type population of sub-L* 
galaxies. From our analysis, we conclude the following. 

• We find no evidence that the slope of the FP 
evolves between Coma and the z = 0.83 cluster 
MS 1054-03. Our sample containing galaxies with 
a ~ 160 km s _1 have a FP with the same slope as 
is found at z = 0. 

• The evolution in M/L is observed to be similar to 
that seen by some previous authors (W04; Holden 



et al. 2005; vv07), where the M/L ratio of the 
cluster population evolves with A log 10 M /Lb = 
-0.50 ± 0.03 from z = 0.023 and z = 0.83, or 
dlog 10 M/L B = -0.60 ± 0.04 dz, with only a mild 
a dependence on the amount of M/L evolution. 
This in marked contrast to more rapid evolution in 
field studies, such as T05 or vdW05, or some other 
cluster results such as J0rgenscn et al. (2006). 



• Using the evolution in the M/Lb and in the col- 
ors of the galaxies, we find a formation epoch of 

= 1.8±oJ to z+ = 2.5tol , the range de- 
pending on the stellar population assumptions we 
make. This is similar to previous measurements, 
see for example iKelson et al. (j2000d ): IWuvts et al.1 
(j200l : iHolden et al.1 (|2005b| ). and vv07 and is in 
goo d agreement with color estimates such as B06 
and lMei et ail (|2009l ). 

• The lack of evolution in the tilt of the FP and in 
the color-cr relations implies little trend in the for- 
mation epoch of cluster galaxies with a. We find 
z* = 2.31q 3 f° r * ne mos t massive, > 300 km s — 
cluster galaxies, while we find = 1.7 jig \ for those 
with ~ 160 km s _1 . This same lack of evolution in 
the tilt of the FP and color-cr relation also implies 
a similar IMF for the lower mass cluster galaxies 
as is seen at higher masses, with a upper limit of 
Ax ~ 0.4 over the mass range we cover, where x is 
the slope of the IMF. 
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APPENDIX 

MS 1054-03 SAMPLE 

Fitting Surface Brightness Profiles 

For each galaxy, we fit the surface brightness profiles in multiple steps. A Sersic model was fit first, using values 
from a catalog derived by SExtractor as initial guesses but with a restriction on n of 1 < n < 4. The output of 
the Sersic model was used an initial guess for the de Vaucouleur's model fit, and the Sersic model was refit using the 
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previous iteration as an initial guess. If the best fitting Sersic index is n — 1, we refit freezing n = 1. This was done 
to ensure convergence. 

Galaxies in crowded environment, generally those near the center of the cluster, were fit simultaneously along with 
all of those galaxies' neighbors. We use the sizes from the SExtractor catalog to find galaxies with close neighbors. 
All galaxies within a radius twice the size of a target galaxy were fit along with the target galaxy. In addition, 
galaxies outside of that radius but with sizes such that the target galaxy lies within twice their size were also included 
in the joint fit. This later criterion was chosen because particularly large galaxies, like the brightest galaxy of the 
cluster, often have extended surface brightness profiles. After the fitting process, if the radii of the target galaxy 
were significantly larger than the initial SExtractor guess, we refit including more neighboring galaxies. In all cases, 
neighboring galaxies were fit with a Sersic model constrained to lie in 1 < n < 4. 

We note here that the ACS magnitudes and surface brightnesses have been adjusted by the offsets as we computed 
from simulations that are listed in Section 12.21 

Errors on the Morphologies 

A number of galaxies in our sample were m orphologically typed m ultiple times. All of the galaxies in our sample 
were classified by Marc Postman twice, once for lPostman et afl (|2005f ) and once for the parent sampl e that the galaxie s 
in this paper are drawn from. In addition, many of these galaxies were classified in WFPC2 imaging (jTran et alj|200l . 
For the purposes of this paper, the main classification is the distinction between a late-type or an early-type system, 
where the former are spiral and irregular galaxies while the later are ellipticals and SO galaxies. The comparison in 
iPostman et al.1 (|2005[ ) finds that 10% of the galaxies in the ACS sample are identified as spirals in the WFPC2 or 
vice- versa. When comparing that two separate classifications done by Marc Postman, we find that 92% of galaxies are 
grouped in the same two broad bins. 

Summary of the Data 

We list all of the galaxies targeted in Table G3 For each galaxy, we list the redshift and dispersion, if available. For 
a summary of why some galaxies do not have dispersions, see Section 12.3.21 The galaxies in the sample of W04 are 
listed with a "w" in front of the identification number. The galaxies with redshifts are those we re-observed as part of 
our program. All of t he other gal a xies h ave spectroscopic information from W04 which does not list redshifts, though 
those can be found in lTran et al.1 ()2007t ). 

Table 3 

Data for Galaxies in MS 1054-03 



ID z log loCT _ log 10 r e a g a n a -OA/mb* B z (U-V) z h r z T J I — J 

(log-LQ km s -1 ) (log 10 kpc) (-0.4 mag) (mag) (mag) (mag) (mag) (mag) 
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838 


2 


329 




0.030 


0.210 


± 


0.008 


0.56 


4.0 


-8.07 


22.78 


1.09 


22.07 


-1 


20.39 


1 


22 


wll92 





840 


2 


156 


: 

± 


0.026 


0.395 


± 


0.005 


0.75 


4.0 


-8.05 


21.80 


1.34 


20.70 


-4 


19.02 


1 


15 


„,i a A n 

wio4y 






2 


386 


■ 


0.050 


0.541 


± 


0.006 


0.78 


4.0 


-8.31 


21.72 


1.34 


20.94 


-4 


19.26 


1 


21 


w2409 






2 


458 




0.050 


0.408 


± 


0.009 


0.68 


3.2 


-8.29 


22.34 


1.43 


21.56 


-4 


19.87 


1 


26 


w3058 




832 


2 


496 




0.028 


1.207 


± 


0.010 


0.72 


3.7 


-9.15 


20.48 


1.31 


19.96 


-5 


18.27 


1 


28 


w3768 






2 


346 


: 

± 


0.047 


0.427 


± 


0.004 


0.90 


4.0 


-8.21 


22.04 


1.26 


21.43 


-5 


19.74 


1 


30 


w3910 






2 


470 


± 


0.062 


0.163 


± 


0.007 


0.42 


4.0 


-7.73 


22.15 


1.49 


21.58 


-2 


19.91 


1 


18 


w4o4o 






2 


526 


: 

± 


0.044 


0.546 


± 


0.009 


0.68 


3 . 4 


-8.28 


21.62 


1.38 


20.74 


-5 


19.05 


1 


28 


w4ozU 






2 


508 


■ 

± 


0.040 


1.339 


± 


0.013 


0.64 


4 . 9 


-9.22 


20.00 


1.40 


19.72 


-5 


18.03 


1 


34 


w4705 






2 


403 




0.062 


0.455 


± 


0.007 


0.28 


1.6 


-8.29 


22.11 


0.73 


21.15 


-1 


19.51 





89 


y\ HOC 

w4y2b 






2 


491 


i 

± 


0.053 


0.267 


± 


0.007 


0.46 


4.0 


-8.00 


22.32 


1.42 


21.48 


-3 


19.78 


1 


35 


w5280 






2 


413 


■ 

± 


0.052 


0.335 


± 


0.011 


0.36 


2.3 


-8.14 


22.31 


1.34 


21.42 


-2 


19.73 


1 


23 


w5298 






2 


453 


± 


0.060 


0.331 


± 


0.012 


0.51 


2.5 


-8.24 


22.60 


1.22 


21.76 


-2 


20.08 


1 


19 


w5347 






2 


405 


± 


0.041 


0.234 


± 


0.008 


0.41 


2.3 


-7.78 


21.93 


1.10 


20.66 


-1 


19.00 


1 


08 


w5450 






2 


369 


± 


0.048 


0.904 


± 


0.005 


0.79 


1.0 


-8.78 


21.09 


1.25 


20.29 


-5 


18.60 


1 


29 


w5529 






2 


260 


± 


0.055 


0.349 


± 


0.009 


0.70 


2.8 


-8.14 


22.25 


1.27 


21.40 


-5 


19.72 


1 


22 


w5577 






2 


484 


± 


0.057 


0.320 


± 


0.016 


0.49 


3.2 


-8.07 


22.21 


1.28 


21.59 


-4 


19.90 


1 


26 


w5666 






2 


456 




0.035 


0.685 


± 


0.005 


0.67 


1.0 


-8.45 


21.34 


1.29 


20.74 


-1 


19.06 


1 


18 


w5756 




831 


2 


379 


± 


0.026 


0.577 


± 


0.004 


0.98 


4.0 


-8.43 


21.83 


1.30 


21.28 


-5 


19.59 


1 


29 


w5840 






2 


405 


± 


0.038 


0.181 


± 


0.007 


0.60 


1.0 


-7.97 


22.66 


1.17 


21.04 


-5 


19.36 


1 


17 


w6036 






2 


396 




0.042 


0.478 


± 


0.012 


0.82 


3.6 


-8.30 


22.01 


1.38 


21.22 


-1 


19.52 


1 


35 


w6301 






2 


438 


± 


0.059 


0.438 


± 


0.004 


0.80 


4.0 


-8.18 


21.91 


1.31 


21.19 


-5 


19.50 


1 


23 



a These are the measured values from fitting a elliptical Sersic model to the two-dimensional images. 

k U — V z refers to the U — V color as if measured by U and V filters redshifted to z = 0.831, the redshift of the cluster. 

c T refers to the morphological class. The scheme is the same as Postman ct al. (2005). E —5 < T < —3 , SO —2 < T < —X, with Spirals 

< T < 10. 



L* Cluster Galaxy M/L Evolution 



21 



B. COMA DATA 

Here we list the compilation of Coma data. These are galaxies in the JFK96 catalog with the corresponding selection 
of E and SO galaxies above a r magni tude limit of r < 15.1 mag AB. Each galaxy has a surface brightness profile derived 
from reprocessed SDSS imaging, see IHolden et all (l200l . This processing mimics what was done for MS 1054-03, so 
we have compara ble data sets. 

As discussed in Hol den et al.l (|2007l ). each galaxy has a B z magnitude derived from the observed g and g — r. Our 
B values match the v alues from JFK96 and w e find good agreement between our B — V colors derived from the g — r 
imaging and those of lEisenhardt et ail (|2007t ). 

The SDSS u imaging is, however, shallow, in comparison to that of lEisenhardt et aD (|2007t ). Therefore, we use the 
colors from that paper. The colors w e list below are e x tincti on corrected and in the redshifted (U — V) 7 , filters. We 
start w ith the colors from Table 8 of Eiscnh ardt et al.l (|2007f ). We use the extinction corrections from Schlcgcl et al. 
(|T998l) . namely A v = 0.052 and A v = 0.032. We redshift the U —V color with the relation (U - V) z = (U - V) - 
0.067(1/ - V) 2 + 0.220([/ - V) - 0.220 as noted in Section^! 

All dispersions come from JFK96. These were measured in equivalent manner as the ones at z — 0.83. 

Table 4 

Data for Galaxies in Coma 



ID log 10 cr log 10 r e q n -0A[i B B z (U - V) z 

(log 10 km s" 1 ) (log 10 kpc) (-0.4 mag) (mag) (mag) 



024 


2.366 


0.290 


0.41 


4.0 


-8.59 


15.03 




027 


2.015 


0.228 


0.70 


2.9 


-8.87 


16.05 




031 


2.416 


1.230 


0.58 


4.0 


-9.57 


12.79 




046 


2.383 


0.519 


0.68 


4.0 


-8.78 


14.36 




049 


2.418 


0.757 


0.54 


4.0 


-8.97 


13.66 




057 


2.232 


0.431 


0.65 


3.4 


-8.71 


14.63 




058 


2.263 


0.805 


0.77 


4.0 


-9.40 


14.50 




065 


2.091 


0.364 


0.78 


2.8 


-9.02 


15.74 


1.48 


067 


2.215 


0.014 


0.91 


3.2 


-8.44 


16.04 


1.44 


068 


2.138 


0.564 


0.95 


4.0 


-9.19 


15.17 


1.46 


069 


2.299 


0.450 


0.82 


3.9 


-8.85 


14.89 


1.60 


070 


2.199 


0.313 


0.96 


4.0 


-8.78 


15.40 


1.48 


072 


2.165 


0.179 


0.48 


2.9 


-8.55 


15.50 


1.39 


078 


2.293 


0.645 


0.83 


4.0 


-9.01 


14.31 




081 


2.183 


0.456 


0.72 


4.0 


-9.16 


15.64 




087 


1.925 


0.065 


0.84 


3.2 


-8.76 


16.59 


1.28 


088 


2.443 


0.093 


0.56 


4.0 


-8.36 


15.44 


1.56 


098 


2.200 


0.332 


0.88 


3.6 


-8.79 


15.33 




101 


2.140 


0.162 


0.57 


3.5 


-8.64 


15.81 


1.43 


103 


2.356 


0.349 


0.53 


4.0 


-8.66 


14.92 


1.50 


104 


2.301 


0.195 


0.37 


3.6 


-8.58 


15.50 


1.51 


105 


2.323 


0.600 


0.69 


3.8 


-8.95 


14.40 


1.58 


106 


2.247 


0.063 


0.89 


3.5 


-8.56 


16.10 


1.32 


107 


1.852 


0.256 


0.82 


2.1 


-8.93 


16.06 


1.30 


109 


2.276 


0.329 


0.95 


3.9 


-8.80 


15.38 


1.58 


116 


2.143 


0.593 


0.49 


4.0 


-9.33 


15.38 


1.46 


118 


2.237 


0.478 


0.87 


3.5 


-8.92 


14.92 


1.51 


119 


2.223 


0.246 


0.56 


4.0 


-8.85 


15.90 


1.37 


120 


2.164 


0.387 


0.82 


3.4 


-8.62 


14.62 


1.46 


121 


2.340 


0.077 


0.70 


2.2 


-8.54 


15.97 


1.47 


122 


1.992 


0.556 


0.86 


4.0 


-9.21 


15.26 


1.29 


124 


2.277 


0.272 


0.47 


2.7 


-8.62 


15.19 


1.48 


125 


2.267 


-0.061 


0.70 


4.0 


-8.33 


16.14 


1.35 


128 


2.065 


0.434 


0.93 


4.0 


-9.17 


15.76 


1.27 


129 


2.376 


1.644 


0.74 


4.0 


-10.01 


11.81 


1.71 


130 


2.352 


0.489 


0.58 


4.0 


-8.88 


14.78 


1.50 


131 


2.251 


0.789 


0.56 


4.0 


-9.33 


14.38 


1.57 


132 


2.142 


0.505 


0.85 


4.0 


-9.36 


15.89 


1.42 


133 


2.375 


0.264 


0.62 


4.0 


-8.55 


15.08 


1.46 


135 


1.920 


0.140 


0.38 


3.1 


-8.92 


16.62 


1.40 


136 


2.221 


-0.172 


0.80 


2.5 


-8.19 


16.34 


1.37 


137 


2.245 


0.364 


0.62 


3.6 


-8.71 


14.96 




113 


2.344 


0.826 


0.73 


3.3 


-9.25 


14.00 


1.68 


111 


2.238 


0.596 


0.49 


4.0 


-9.07 


14.70 


1.48 


115 


2.152 


0.415 


0.39 


3.3 


-9.01 


15.47 


1.48 


146 


2.042 


0.549 


0.63 


3.2 


-9.24 


15.37 


1.37 


148 


2.572 


1.234 


0.39 


4.0 


-9.33 


12.18 


1.75 


150 


2.044 


0.353 


0.42 


4.0 


-8.94 


15.61 


1.40 


151 


2.184 


0.784 


0.66 


4.0 


-9.32 


14.39 


1.42 


152 


2.217 


0.616 


0.63 


4.0 


-9.27 


15.10 


1.52 


153 


2.153 


0.359 


0.89 


4.0 


-8.97 


15.65 


1.44 


155 


2.203 


0.783 


0.74 


4.0 


-9.35 


14.48 


1.50 


156 


2.025 


0.268 


0.59 


3.6 


-9.08 


16.36 


1.29 


157 


2.142 


0.305 


0.54 


4.0 


-8.97 


15.92 


1.32 
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Data for Galaxies in Coma 



159 


2.306 


0.556 


0.37 


3.4 


-8.88 


14.43 


1.54 


160 


2.284 


0.513 


0.68 


3.9 


-9.04 


15.05 


1.50 


161 


2.251 


0.547 


0.53 


4.0 


-8.95 


14.65 




167 


2.336 


0.509 


0.74 


3.9 


-8.81 


14.49 


1.53 


168 


2.347 


0.264 


0.93 


2.5 


-8.48 


14.88 


1.53 


170 


2. 174 


0.445 


0.52 


3.7 


-9.00 


15.30 


1.44 


172 


2.227 


0.143 


0.63 


3.6 


-8.52 


15.60 


1.45 


173 


2.179 


0.161 


0.88 


a n 
4.0 


-8.66 


15.87 


1.41 


1 T A 

174 


2.28 / 


0.014 


U.24 


2.3 


-8.28 


15.64 


1.50 


1 

1 / 5 


2.263 


0.378 


0.62 


2.9 


-8.78 


15.08 


1.58 


176 


2.249 


0.146 


0.87 


1.6 


-8.60 


15.78 


1.50 


177 


2.038 


0.218 


0.18 


3.7 


-8.93 


16.25 


1.40 


1 /9 


O Af\A 

Z.4U4 


0.478 


0. /4 


a n 
4.0 


-8.63 


1 A OA 

14. zU 


1.51 


181 


2.194 


0.089 


0.96 


3.9 


-8.60 


16.07 


1.38 


191 


1.994 


0.007 


0.87 


2.3 


-8.54 


16.32 




192 


1.989 


0.178 


0.23 


1.6 


-8.65 


15.75 




193 


2.101 


0.155 


0.38 


2.9 


-8.74 


16.08 




194 


2.398 


0.646 


0.68 


4.0 


-8.93 


14.10 




204 


2.116 


0.416 


0.31 


4.0 


-9.06 


15.57 




206 


2.351 


0.627 


0.28 


3.3 


-8.83 


13.96 




207 


2.192 


0.181 


0.24 


2.6 


-8.67 


15.79 




210 


2.246 


0.153 


0.51 


3.4 


-8.53 


15.58 




217 


2.317 


0.766 


0.81 


4.0 


-9.18 


14.12 




238 


2.041 


0.037 


0.86 


3.4 


-8.51 


16.10 




239 


2.359 


0.791 


0.62 


4.0 


-9.18 


14.01 




240 


2.415 


1.024 


0.61 


4.0 


-9.31 


13.18 





